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The associated production of neutral Higgs bosons with the Z gauge boson (h°Z°, H°Z°) is in- 
vestigated in the context of the future linear colliders, such as the ILC and CLIC, within the general 
two-Higgs-doublet model (2HDM). We compute the corresponding cross-section for the processes 
e + e -^Z h /Z°H at one-loop, including the full set of corrections at 0(a^ w ) together with the 
leading 0(a^ m ) terms, in full consistency with the available theoretical and phenomenological con- 
straints. We find that the wave- function corrections to the external Higgs fields are the dominant 
source of the quantum effects, which turn out to be large and negative (e.g. S a jo ~ —20% / — 60%) 
in all the y/s range, and located predominantly in the region around tan f3 ~ 1 and moderate values 
of the parameter A5 (being A5 < 0). This behavior can be ultimately traced back to the enhance- 
ment potential of the triple Higgs boson self-couplings, a trademark feature of the 2HDM with no 
counterpart in the Higgs sector of the Minimal Supersymmetric Standard Model. Even under this 
substantial depletion of the one-loop corrected signal (which is also highly distinctive with respect 
to the SM expectation for e + e~ — > Z°H), the predicted Higgs-strahlung rates comfortably reach a 
few tens of femtobarn, which means barely ~ 10 3 — 10 4 events per 500 fb" 1 of integrated luminosity. 
Due to their great complementarity, we argue that the combined analysis of the Higgs-strahlung 
events (h°Z°, H°Z°) and the previously computed one-loop Higgs-pair production processes (h°A°, 
H°A°) could be instrumental to probe the structure of the Higgs sector at future linac facilities. 

PACS numbers: 12.15.-x, 12.60.Fr, 12.15.Lk 



I. INTRODUCTION 

After more than 40 years since the seminal ideas 
coined by a handful of theoretical pioneers [l| , our present 
understanding of the Electroweak Symmetry Breaking 
(EWSB) phenomenon through the Higgs (Englert-Brout 
and Guralnik-Hagen-Kibble) mechanism is still rather in- 
complete and experimentally inconclusive. On the one 
hand, there is no compelling alternative to consistently 
embed the EWSB mechanism into the quantum field the- 
oretical description of particle physics offered by the - in 
so many regards successful - Standard Model (SM) of 
the Strong and Electroweak interactions. On the other 
hand, we do not have a single phenomenological hint of 
the existence of elementary scalar fields, not to mention 
the fact that we do not understand how to make com- 
patible the EWSB mechanism and its associated vacuum 
energy with fundamental problems of different scope, for 
example in the domain of cosmology. Still, the possibil- 
ity to describe the inner theoretical structure of the SM 
through the EWSB mechanism is so successful within 
the restricted particle physics domain that it would not 
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be wise, not even advisable, to cease our pursue of the 
phenomenological implications of the Higgs mechanism 
paradigm till the limits of the current experimental pos- 
sibilities. 

It goes without saying that the quest for experimen- 
tal evidences of the Higgs boson is a most preeminent 
milestone of the upcoming generation of collider facili- 
ties. Nonetheless, the future data might well reveal that 
the purportedly found Higgs boson actually belongs to a 
richer model structure, which might be grounded some- 
where beyond the minimal conception of the SM, namely 
of a single, spinless, fundamental constituent of matter. 
If so, a few more fundamental spinless constituents could 
appear. A particularly well-motivated extension is the 
two-Higgs-doublet model structure encompassed by the 
Minimal Supersymmetric Standard Model (MSSM) 0. 
Here the physical Higgs boson spectrum contains a couple 
of CP-even (h°, H°), one CP-odd (A ) and two charged 
Higgs bosons (H ± ), and the corresponding Higgs poten- 
tial is highly constrained by the underlying Supersym- 
metry (SUSY). One particular consequence of the latter 
is that SUSY invariancc restricts all Higgs boson self- 
interactions to be gauge-like. While this represents an 
economy from the point of view of the number of cou- 
plings in the potential, it is strongly imbalanced by the 
exceeding number of parameters (mixings and masses) 
populating the other sectors of the MSSM. In addition, 
the pure gauge nature of the Higgs boson self-couplings 
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make them highly inconspicuous from the practical point 
of view, in the sense that they are unable to trigger any 
outstanding phenomenological signature. The core of the 
enhancement capabilities of the MSSM Lagrangian re- 
sides, instead, in the multifarious pattern of Yukawa cou- 
plings between the Higgs bosons and the quarks, as well 
as between quarks, squarks and charginos/neutralinos. 
The rich interplay of opportunities that they give rise 
to has been extensively analyzed in the past within a 
plethora of processes, see e.g. 0-0], and also [Hill for 
reviews on the subject. 

The two-Higgs S , L r i(2)-doublet structure of the Higgs 
sector in the MSSM constitutes a genuine prediction of 
the SUSY dynamics. Nonetheless, a more general archi- 
tecture can appear in the form of a non-SUSY frame- 
work through the so-called general (unconstrained) two- 
Higgs-doublet Model (2HDM) 1 . Although they share a 
common Higgs boson spectrum, the potentially most dis- 
tinctive phenomenological features of both models are 
located in very different sectors. The most relevant ob- 
servation here is that, in the absence of an underlying 
SUSY, the triple (3H) and quartic (4H) Higgs-boson self- 
interactions are no longer restrained to be purely gauge. 
This can have a tremendous impact e.g. in the physics 
of the top quark in hadron colliders, as it was shown 
long ago in [T^], and it can also trig ger significant neu- 
tral flavor-changing interactions 11311 that can perfectly 
compete with the SUSY ones fLU Il5| . 

Much attention has been devoted to Higgs boson pro- 
duction and decay in hadron colliders, extending from 
the ongoing Tevatron facility at Fermilab to the brand 
new LHC collider recently operating at CERN [l-Qjj]. 
However, precision Higgs boson physics will greatly ben- 
efit from the interplay |l8| of the forthcoming generation 
of linear colliders (linac), such as the ILC and CLIC 
projects [IH . Here a variety of processes can provide new 
clues to Higgs boson physics, e.g. the production of triple 
Higgs-boson final states, both in the MSSM [2(| and in 
the general 2HDM 2l| ; the double Higgs-strahlung chan- 
nels hhZ° [22j]; and the inclusive Higgs- pair production 
via gauge-boson fusion [2(|[23j]. In the same vein, also the 



77 mode of a linac has been explored [24j , in particular 
the loop- induced production of a single neutral Higgs bo- 
son [25[ and of a Higgs boson pair [2g|. In all the above 
mentioned cases, promising signatures were singled out 
and illustrate that, if effectively realized in Nature, such 
hints of a generic 2HDM structure could hardly be missed 
in the superbly clean environment of a linac. And, what 
is more, they could not be confused as having a SUSY 



1 To be more precise, we refer to this model as "general" in the 
sense that we allow all possible operators leading to a rcnormal- 
izablc, gauge-invariant and CP-conserving Higgs potential. As 
usually done in practice , we impose an additional (softly- 
broken) Z2 discrete symmetry <$>i — ¥ (— 1)*$^ where i = 1,2 
denote each of the SUl{2) Higgs doublets, as a sufficient condi- 
tion to banish the tree-level flavor changing neutral currents. 



origin, because of the intrinsically different nature of the 
Higgs self-interaction sector. Besides, outstanding finger- 
prints of a generic 2HDM could also be stamped in the 
pattern of radiative corrections to Higgs production pro- 
cesses; for instance, quantum effects on the cross-sections 
for two-body Higgs boson final states: 

e+e~^2h (2/i = li°A°; H°A°; H+H~), (1) 



have been attentively investigated in the MSSM [27H30| . 
As for the general 2HDM, the efforts were first concen- 
trated on the production of charged Higgs pairs [3lJ . This 
program has been recently brought to completion from 
a full-fledged study of the quantum effects on the pro- 
duction cross-sections in the neutral Higgs sector [321 ] . 
Alongside these processes we find the more traditional 
Higgs-strahlung events, in which a Higgs boson is pro- 
duced in association with the Z°: 



e ' e 



h = h°,E°. 



(2) 



(Notice that the A°Z° final state is forbidden by CP- 
conservation.) Processes Q are complementary to the 
(TT|) ones - cf . H3 , and also [2(], H3 for a phenomeno- 
logical analysis in the MSSM context. Let us also re- 
call that the last available limit on the SM Higgs boson 
mass, placed by LEP searches, comes precisely from in- 
vestigating the "Bjorken process" [35| , i.e. e + e~ — > H Z° 
followed by Z° //, with the result: M H > 114.4 
GeV [3(|. In this article, we aim at discussing the corre- 
sponding one-loop corrections to the generalized Bjorken 
processes within the 2HDM. However, in contrast to 
the original SM case, where the radiative corrections are 
small, the quantum effects on the processes © can be 
large and are mainly driven by the 3H self-couplings. In 
fact, our main aim here is to identify the regions of the 
parameter space where this is so, and then quantify the 
impact of the potentially enhanced 3H self-couplings on 
the final cross-sections. In this way, joining this study 
with that of Ref. [HJ on the pairwise Higgs boson pro- 
duction channels (Q} in the 2HDM, a rather complete 
panorama of the genuine quantum effects associated to 
the basic neutral Higgs boson production channels be- 
comes available. 



II. HIGGS-STRAHLUNG EVENTS AT 
ONE-LOOP: THEORETICAL SETUP 

A good deal of attention has been devoted in the lite- 
rature to multiple properties of the Higgs-strahlung pro- 
cess e + e~ — > HZ in the SM, see e.g. [37| and references 
therein. Not surprisingly it was one of the gold-plated 
channels for Higgs boson search at LEP. Concerning the 
MSSM, the Higgs-strahlung channel s dp have also been 
discussed extensively - cf. Refs. [H, [381] and part two of 
the review Q - and are currently under investigation also 
in the MSSM with CP-violating phases [3!| (the so-called 
complex MSSM [4C|). The upshot of these analyses spot- 
lights the following features: i) the dominant source of 
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corrections at one-loop originate from the Higgs boson 
propagators, which can be reabsorbed into an effective 
(loop-corrected) mixing angle a e ff (equivalently, a more 
generic mixing matrix in the complex MSSM case); ii) 
the corrections to the ZZ h vertex are in general small, 
although they can reach the level of 10 % for very low 
(or high) values of tan /3, precisely in the regions where 
the Higgs Yukawa couplings to heavy quarks become en- 
hanced; iii) the electromagnetic corrections to the initial 
state with virtual photonic corrections and initial-state 
radiation do not differ from the SM case, being in general 
large and positive (except near the production thresh- 
old). In the present article, our endeavor is to seek for 
the genuine phenomenological imprints associated to the 
generic (unconstrained) 2HDM dynamics, most particu- 
larly to the potentially enhanced 3H self-couplings. 

Let us recall that the general 2HDM [ll[ is obtained 
by canonically extending the SM Higgs sector with a sec- 
ond SUl(2) doublet with weak hypercharge Y = +1, so 
that it contains 4 complex scalar fields. The free param- 
eters Xi in the general, CP-conserving, 2HDM potential 
can be finally expressed in terms of the masses of the 
physical Higgs particles (M^o, Mjj°, M^o, Mff±), tan/3 
(the ratio of the two VEV's (H®) giving masses to the 
up- and down- like quarks), the mixing angle a between 
the two CP-even states and, last but not least, the self- 
coupling A5, which cannot be absorbed in the previous 
quantities. Therefore we end up with a 7-free param- 
eter set, to wit: (M h o, M H o, M^o, M H ±, sin a, tan/3, 
A5). Furthermore, to ensure the absence of tree- level fla- 
vor changing neutral currents (FCNC), two main 2HDM 
scenarios arise: 1) type-I 2HDM, in which one Higgs dou- 
blet couples to all quarks, whereas the other doublet does 
not couple to them at all; 2) type-II 2HDM, where one 
doublet couples only to down-like quarks and the other 
doublet to up-like quarks. The MSSM Higgs sector is 
actually a type-II one, but of a very restricted sort (en- 
forced by SUSY invariance) [ll[ . We refer the reader to 
Ref. [32[ for a comprehensive account on the structure 
of the model and for notational details. In particular, 
Table II of that reference includes the full list of trilinear 
couplings within the general 2HDM that are relevant for 
the present calculation. Further constraints must be im- 
posed to assess that the SM behavior is sufficiently well 
reproduced up to the energies explored so far. Most par- 
ticularly, we take into account i) the approximate SU (2) 
custodial symmetry, which can be reshuffled into the con- 
dition \Sp2HDAi\ < 10~ 3 [4l|; and ii) the agreement with 



the low-energy radiative S-meson decay data (which de- 
mand M H ± > 300 GeV for tan/3 > 1 HH in the case of 
typc-II 2HDM). 

Additional requirements ensue from the theoretical 
consistency of the model, to wit i) perturbativity [43l[: ii ) 
unitarity [44j ; and iii) stability of the 2HDM vacuum [45| ■ 
Let us expand a bit more on the latter conditions, which 
turn out to play a crucial role in our analysis. As for 
the unitarity bounds, and following Ref. [4J|: the basic 
underlying strategy is to compute the S-matrix elements 



Sij for the possible 2 — s- 2 processes involving Higgs and 
Goldstone bosons in the 2HDM, and to subsequently re- 
strain the corresponding eigenvalues {/,/- SuUf- 1 = a.; 5^ 
by the generic condition | on \ < 1/2 Vi. The latter require- 
ment translates into a set of upper bounds on the quartic 
scalar couplings Aj, i = 1 . . . 6 - and hence on combina- 
tions of Higgs masses, trigonometric couplings and, most 
significantly, the parameter A5. As far as vacuum stabil- 
ity restrictions are concerned [45j, they may be written 
as follows: 



Ai + A 3 > ; A 2 + A 3 > 
2 v / (Ai+A 3 )(A 2 +A3) + 2A 3 + A 4 
+min[0, A5 — A4, A6 — A4] > 0. 



(3) 



More refined versions of these equations may be obtained 
from the Renormalization Group (RG) running of the 
quartic couplings \i(fJ?) at high energies. Nonetheless, 
for our current purposes we need not assume any par- 
ticular UV completion of the 2HDM. This would be an 
unnecessary additional assumption at this stage of the 
phenomenological analysis of our processes. Therefore, 
we are not tied to any specific UV-cutoff (which could 
be, in principle, as low as A ~ 1 — lOTcV). In this 
sense, we may just apply Eq. <j3j> with all Ai taken at the 
EW scale. This is after all the scale at which we fix all 
our input parameters and perform the phenomenologi- 
cal analysis (including the renormalization) of the pro- 
cesses under study. By the same token, the (cutoff de- 
pendent) triviality bounds are also circumvented. The 
latter kind of restrictions only impose tight upper limits 
on the Higgs boson masses when a very large UV-cutoff 
(e.g. the Planck scale or the GUT scale) is considered. 
Our choices of Higgs boson mass spectra (cf. Tabic H]) 
lie, in any case, in a sufficiently low mass range so as to 
conform with the typical mass requirements allowed by 
triviality J4(| - even for large cutoff scales. More details 
on the constraints setup are provided in Ref. [32j . 

The dynamics of the Higgs-strahlung processes under 
consideration is driven at leading order by the tree-level 
interaction Lagrangians: 

e sin(/3 - a) M z 

£z°z°h° = : — : 9 L n L „ h > 



Sw cw 



£ z o z o H n — 



eoofl^-a)Jlfe_ irz o z o H | (4) 



sw Cw 



where sw = sm6*vj/ and c\y = cos 9w for the elec- 
troweak mixing angle 9w- Since £A°h°z ^ cos(/3 — a) 
and £a°h°z° °t sin(/3 — a), it is clear that the processes 
([2|) arc complementary to the Higgs pair production ones 
((T|). Furthermore, as all these couplings (U]) are generated 
by the gauged kinetic terms of the Higgs doublets (cf. 
Eq. (30) of Ref. [HI), they are fully determined by the 
gauge symmetry and hence show no intrinsic difference in 
the 2HDM as compared to the MSSM. At the end of the 
day, this is the reason why a tree-level analysis of these 



FIG. 1: Set of Feynman diagrams contributing to e e~ — > h Z at one-loop level within the 2HDM. The shaded blobs stand 
for all possible loops with virtual particles. An equivalent collection of diagrams accounts for the complementary process 
e + e~ — > H°Z°. The last couple of rows displays the complete set of one-loop counterterm diagrams. 



events is most likely insufficient to disclose their true na- 
ture. Similarly, in the limit a = f3 - tt/2 the h°Z°Z° 
coupling coincides with the analogue coupling in the SM, 
HZ°Z°. It is, therefore, the pattern of radiative correc- 
tions the characteristic signature associated to each one 
of the possible models; most particularly, it should help 
to disentangle SUSY versus non-SUSY extended Higgs 
physics scenarios. 

The leading-order 0(a ew ) scattering amplitude follows 
from the s-channel Z°-boson exchange and renders 



7U(o)(e+ e -->h Z°) = - 



e 2 M z sin(/3 - a) 



sw cw (s — M§) 

X^(Pl,m) f- (&2,CT 2 ) (gL PL + 9RPR)u(p 2 ,m) 



(5) 



Here and 7712 refer to the 4-momenta and helicitics 
of the electron and positron, and e (fe, 02) is the polar- 
ization 4- vector of the Z gauge boson with 4-momentum 
&2 and helicity a%. We have introduced also the left and 
right-handed weak couplings of the Z° boson to the elec- 
tron, g L = (-1/2 + s 2 w )/c w s w , g R = s w /c w , and the 
left and right-handed projectors Pl.r = (1/2)(1 =F 75). 
Let us notice that we do not include the finite Z°-width 
corrections, since they are completely negligible for the 
ccnter-of-mass energies that we consider here. Finally, 



the total cross section er(c + e~ — > h°Z°) at the tree-level 
is obtained after squaring the matrix element (O, per- 
forming an averaged sum over the polarizations of the 
colliding e + e~ beams and the outflowing Z° boson, and 
integrating over the scattering angle. 

The calculation of cr(e + e~ — > h°Z°) at one loop is 
certainly much more cumbersome. To start with, it is 
UV-divergent and it thus requires of a careful renormal- 
ization procedure in order to render finite results. We 
adopt here the conventional on-shcll scheme in the Feyn- 
man gauge [13]. In the MSSM, these cross-section cal- 
culations can be mostly carried out in an automatized 
fashion through standard algebraic packages which allow 
a rapid and efficient analysis of the electroweak preci- 
sion observables (48j . Several public codes are available, 
see e.g. (49l - l52j . However, in our case the calculation is 
non-supersymmetric and we must deal with the renor- 
malization of the Higgs sector in the class of generic 
2HDM models. A detailed description of the renormal- 
ization procedure for the 2HDM Higgs sector in the on- 
shell scheme has been presented in [32( and we refer the 
reader to this reference for all the necessary details. We 
will make constant use of the framework described in this 
reference and sometimes we will refer to particular for- 
mulae of it. 

With this renormalization setup in mind, the vari- 
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FIG. 2: Set of Feynman diagrams contributing to e + e~ — > 
h°Z° at one-loop level within the 2HDM. These diagrams des- 
cribe the box-type quantum corrections. An equivalent col- 
lection of diagrams accounts for the complementary process 





FIG. 3: Sample of Feynman diagrams describing the one- 
loop corrections to the h°Z°Z° interaction for the process 
e + e~ — > h Z within the 2HDM, including quark-mediated 
(top line) and Higgs boson-mediated quantum corrections 
(bottom line). 



ous contributions to the scattering amplitude of e + e — > 
h°Z° at one-loop (Figs. EEJ can be classified in a mean- 
ingful way. First of all, let us note that all of the one- 
loop diagrams are of course of order 0(a^ w ). However, 
some of them include, in addition, enhancement factors 
sourced by the trilinear Higgs boson couplings A3H, see 
Table II of " 
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In such cases, we shall include these 
factors when assessing the order of magnitude of the di- 
agram. We are now ready for sorting out the one-loop 
contributions in different categories: 

• Self-energy corrections to the Z° and 7 — Z° mixing 
propagators, all of them of order 0(oq w ) with no 
enhancement factors at this order. 

• Vertex corrections to the e + e~Z° interaction, which 
are also of order O(<o%, w ). It should be noted that 
we do not include the virtual photonic 0(a ern a ew ) 
effects nor the real brcmsstrahlung emission off 
the e + /e~ legs. These pure QED corrections and 
the weak ones factorize into two subsets which are 



separately UV finite and gauge- invariant. More- 
over, these photonic contributions are fully insen- 
sitive, at the order under consideration, to the rel- 
evant 3H self-couplings on which we focus. For 
this kind of processes involving electrically neutral 
Higgs bosons in the final state, the one-loop QED 
effects arc confined to the initial e + e~ vertex. In 
practice, the net outcome of the accompanying ini- 
tial state radiation (ISR) is to lower the effective 
center-of-mass energy available for the annihilation 
process. All in all, these effects are unessential at 
this stage to test the presence of the new dynam- 
ical features triggered by the 2HDM in the Higgs- 
strahlung events under analysis. 

• Vertex corrections to the e + e~H interaction. Since 
we have explicitly set m e = throughout our calcu- 
lation, the e + e~H tree-level Yukawa coupling is ab- 
sent and the corresponding one-loop diagrams au- 
tomatically render a UV finite contribution which 
is, in any case, very small. 

• The loop-induced 7Z h° interaction. This one is 
order 0(ea ew X^h) and therefore includes an en- 
hancement factor A3H . Due to the 7 — Z° mixing 
at one-loop, the following counterterm is needed so 
as to render a UV-finite vertex: 
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h n Z° 
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(6) 



and similarly for the counterterm associated to the 
effective £h°z°7 interaction: 



8 £h°z°~ 



e cos(/3 - a) M z 
2swcw 



5Z 7Z og^ZlA v B° 



(7) 



• The vertex correction for hZ°Z° encompasses dif- 
ferent contributions of order 0{cq w , a ew eX^n) 
(see the sample diagrams in Fig. And, most 
important, there are, in addition, corrections of 
0(a ew X^ H ) (in fact, the dominant ones) which are 
characterized by the highly conspicuous enhance- 
ment factor squared X^ H . The latter originates 
from the associated vertex counterterm, most par- 
ticularly from the Higgs field rcnormalization con- 



stant Z, 



1/2 
h n 



1 + (1/2)5 Z\ x a. As we shall discuss in 



more detail below, 6 Z h o is sensitive to the scalar- 
scalar self-energy and hence it involves products 
of two triple Higgs self-couplings (see Fig. [4]). Of 
course, an equivalent discussion holds for the com- 
plementary channel, H°Z°. 

The full form of the associated one-loop vertex 
counterterms ensues from the usual splitting of bare 
fields and parameters into the renormalized ones 
and associated counterterms, to wit: 



£ (.90 = 9 



sc. 



(8) 
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For the h°Z°Z° and H°Z°Z° case one gets respec- 



tively: 
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where we can spot in the structure of these formu- 
lae the presence of the above mentioned Higgs field 
renormalization counterterms S Z^o and S Z^o . The 
term 6 Zh°H° accounts for the WF mixing h° o H° 
emerging from the Z°Z°h° and Z°Z°H° vertices. 

The finite wave-function (WF) correction to the ex- 
ternal Higgs fields h° or H° (Fig.iJ. These are re- 
lated to the fact that we have chosen the residue of 
the A°-propagator at its pole to be one, and there- 
fore there is no more freedom to make the same 
choice for the other Higgs bosons. This entails a 
finite WF renormalization correction, see Ref. 32 1 
for details. These finite renormalization effects are 
of utmost importance in this case, as they trigger 
leading contributions of order 0(a ew \\ H ). They 
actually drive the very bulk of the quantum effects 
(see below for a more detailed discussion). 



• Finally, the box- type 0{a 2 ew ) diagrams (Figure [2]), 
whose contribution is non-negligible, in particular 
at large center-of-mass energies - since they are not 
suppressed as 1/s. 

The dominance of the WF corrections is a very char- 
acteristic feature of the processes under analysis. It orig- 
inates from a non-trivial cancelation between the Higgs 
boson field counterterms, which appear in two different 
pieces of the overall one-loop amplitude, though with op- 
posite signs. Let us further elaborate on this important 
point. Without loss of generality, we concentrate on the 
h°Z° channel. The complete set of contributions to the 
scattering amplitude up to the one-loop level may be split 
in the following manner: 



M 



(o+i) 



>h°Z° 
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c+c-^h°Z° 



► h°Z° 
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wherein the finite WF corrections to the external Higgs 
field arc introduced as 



•^e+e-^hfZ 



>h°Z° 



(0) 

C + C " 



(12) 



At this point we are making explicit use of the on-shell 
renormalization conditions defined in [321 ] . in particular 
3fe ShOH°(A£? ) = 0- Moreover, we only retain those 
contributions that are leading-order in the triple Higgs 
self-couplings. Let us notice that the counterterm piece 
SM^ is sensitive to 0(\\ H ) effects through the WF 
renormalization of the Higgs fields. More specifically, 
from the explicit expression of the h°-field counterterm 
as a linear combination of the two Higgs doublet coun- 
terterms (see Section IV of Ref. [Hf for details) 

5 Z h o = sin 2 a 5Z<s> 1 + cos 2 a 5Z<j> 2 = — 5fte [M\ ) + • • • 

we may single out the following C(A 2 H ) contribution 
from Eq. ©, 



xr _^ e sin(/3 - a) 

L Sw Cw 

e sin(/3 — a) 



2 sw cw s. 



5ie Z' A o{Ml ) + 



0(A| H ) 



(14) 



But, as warned, the finite WF renormalization factor of 



(13) 
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the h° field in (fT2f is sensitive to 0(X 2 H ) terms too: 



5Re £' hQ (M 2 ) = 3?e £ h0 (M 2 „) + S Z 



Ike Y,' h0 (M h 2 ) - 3?e £^ (Af£ ) 



0(A 



(15) 



where the dots stand for those terms with dependencies 
other than 0(X 2 H ). Notice, therefore, that part of the 
0(X 2 H ) dependence cancels out between the counterterm 
diagram associated to the h°Z°Z° vertex - cf. Eq. ((H)) 
- and the finite h° WF- factor - cf. Eqs. (JT2J) and (|T5i 
Specifically, it is the piece 5fte (M^ ) that exactly can- 
cels between the two terms. The remaining 0(X 2 H ) con- 
tributions come from 5ReEj i0 (M^ ) and are generated by 
the Feynman diagrams displayed in Fig. 2J When suffi- 
ciently enhanced, these pieces account for the bulk of the 
one-loop quantum corrections. A typical diagram in the 
class of 2-point functions provides the following contri- 
bution 2 : 



M 



in 



1 



+ e -^h°Z° - ~2 ffieE h°( M li°)^e+e-->.h°Z° 



-lM^ h o z o\M H \ 2 X 



,4-D 



dp 2 v "'./ (2tt)^ (g 2 -A/ 2 )[(g + p) 2 -M 2 ] 
■5 Toff ^Si-^ozo Refi^o, Af 2 , M 2 ) , 

(16) 



wherein the scalar two-point function is defined as in 
Ref. M : 
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/' 



4-D 



(2tt) d (q 2 -mj) [(p + q) 2 -m 2 
B (p 2 ,mj,ml), 



16 7P 



(17) 



/i denoting the 't Hooft mass unit. In the expression 
(fT6"|) . M denotes the typical mass scale(s) which appear 
in these one-loop 2-point functions. Owing to the deriva- 
tive with respect to the external momentum, the 1-point 
functions do not contribute to the WF renormalization, 
and hence the quartic Higgs boson self-couplings are not 
involved in this calculation. This means that the first 
line of diagrams in Fig. 2] does not actually contribute. 
We emphasize that the overall sign for this expression de- 
pends on the sign of such 2-point function. Notice that 



2 Notice that we employ the convention that iS equals the Higgs 
boson self-energy diagram, such that S does not contain the 
global imaginary part emerging from the loop integral. This is 
why we have multiplied the loop integral in l|16|l by —i before 
taking the real part. 



VleB (p 2 ,M 2 ,M 2 ) > for p 2 < 4M 2 , which is the case 
we wish to focus on for the h°Z° channel, and also for 
H°Z° (as long as the resonant decay H° — > h°h° is forbid- 
den by kinematics). Therefore, in most of the scenarios 
of interest, such finite WF corrections carry an overall 
minus sign. Being proportional to X 2 H , they become the 
leading quantum effects in the region where the trilincar 
couplings are enhanced, and as a result they generate 
a characteristic pattern of quantum effects in which a 
systematic suppression of the tree-level cross-section is 
predicted. 

The presence of the large negative corrections induced 
by the Higgs boson self-energies brings forward a charac- 
teristic signature for the production cross-sections of the 
Higgs- strahlung processes $Z§. This feature is in marked 
contradistinction to the situation with the Higgs boson 
pair production mechanisms ([1]), where the correspond- 
ing corrections are just opposite in sign, i.e. large and 
positive, see (32[. We shall further comment on these 
interesting and correlated features in the next section. 
From the general structure of Eq. (fT5|) , one may antici- 
pate the typical (maximum) size of the quantum effects 
on the Higgs-strahlung processes as follows: 



5 r 



a 



-(o+i) _ o-flO _ (2M^M^ 
' (\M^\ 2 ) 



l A 3g| 2 f( M 2 M 2 



! ,M 2 ) 



(18) 



/ being a dimensionless form factor (basically account- 
ing for the behavior of the B' Q function). The notation 
(. . .) stands for the various operations of averaging and 
integration of the squared matrix elements. Taking into 
account that unitarity limits let the trilincar couplings 
reach values as large as jj 11 - ~ ~ 30, and assuming 
M ~ 200 GeV and / ~ 0(1), the above estimate typica- 
lly predicts a strong depletion of the tree-level signal by 
5 r ~ —90%. This prediction falls in the right ballpark 
of the exact numerical results that will be reported in 
the next section, which point to a maximum depletion of 
S r ~ -60%. 

Before closing this section, let us recall that the coun- 
terterm amplitude 5 M. + } _ _vh°z derives from a number 
of renormalization conditions that determine the renor- 
malized coupling constants and fields in a given renor- 
malization framework. As we have said, in our calcula- 
tion the renormalization is performed in the conventional 
on-shell scheme in the Feynman gauge, appropriately ex- 
tended to include the 2HDM Higgs sector. For the latter, 
we need in particular a renormalization condition for the 
parameter tan (3 = V2/V1. We adopt the following [53| : 

5 vi S V2 
Vi v 2 



(19) 



This condition insures that the ratio vi / v\ is always ex- 
pressed in terms of the true vacua after the renormaliza- 
tion of the Higgs potential. The corresponding countert- 
erm resulting from tan j3 — > tan f3 + 5 tan j3 can then be 
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computed explicitly: 
Stan/3 



tan p M z sin 2f3 



5ReS A o z o(Mlo). 



(20) 



This counterterm is involved in Eqs. (|9))- (|T0| . and it 
also determines the WF mixing term S Z h a H a that 
appears in these equations, as follows: S Z h o H o = 
sin 2a (S tan f$/ tan/3). We refer the reader once more to 
the exhaustive presentation of Ref. [32| for the renormal- 
ization details. 
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due to their mass splittings (and also to the fact that 
M h o > 140 GeV), Sets C and D can only be realized in 
a general (non-SUSY) 2HDM framework. Furthermore, 
in view of the value of the charged Higgs mass, Sets A-C 
are only suitable for type-I 2HDM, whereas Set D is valid 
for cither type-I and type-II 2HDM's. 

Apart from reflecting a variety of possible situations 
in the 2HDM parameter space, some of these sets can be 
mimicked by the Higgs boson mass spectrum in super- 
symmetric theories. For instance, Sets A and B can be 
ascribed to characteristic benchmark scenarios of Higgs 
boson mass spectra within the MSSM; in particular, Set 
B lies in the class of the so-called maximal mixing scenar- 
ios [54[, for which AI h o takes the highest possible values 
within the MSSM. The numerical mass values for the 
MSSM-like sets have been obtained with the aid of the 
program FeynHiggs by taking the full set of EW correc- 
tions at one- loop [50[ • 

We remark that the more massive the Higgs bosons 
are, the stronger the constraints that unitarity imposes 
over | A5 1 . The maximum (negative) values are roughly 
attained for A 5 ~ -9 (Set A); A 5 ~ -10 (Sets B and C); 
and A 5 ~ -8 (Set D). 
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FIG. 4: Subset of Feynman diagrams driving the dominant 
contribution to the h° boson self-energy. Notice that the two- 
point functions are sensitive to two triple Higgs-boson self- 
couplings. The tadpole diagrams, however, which depend on 
quartic couplings, do not contribute to the WF renormaliza- 
tion. An equivalent collection of diagrams contribute to the 
corresponding H° self-energy. 



III. NUMERICAL ANALYSIS 

In this section, we present the numerical analysis of the 
one-loop computation of the Higgs-strahlung processes 
([2]) within the 2HDM. We shall be concerned basically 
with the following two quantities: i) the predicted cross- 
section at the Born-level a^ 1 and at one loop ct^ 0+1 ^; and 
ii) the relative size of the one-loop quantum corrections: 



fir 



r (0+l) _ a (o) 



(21) 



We have carried out our analysis with the help of the 
standard algebraic and numerical tools FeynArts. Form- 
Calc and LoopTools (49| . 

The different Higgs boson mass sets that shall be used 
hereafter are quoted in Table HI Let us highlight that, 





M h o [GeV] 


M H o [GeV] 


M A o [GeV] 


M H ± [GeV] 


Set A 


115 


220 


220 


235 


Set B 


130 


160 


150 


170 


Set C 


140 


150 


200 


200 


Set D 


150 


200 


260 


300 



TABLE I: Choices of Higgs masses (in GeV) that are used 
throughout the calculation. Due to the values of the M H ± 
mass, Sets A-C would only be suitable for type-I 2HDM, while 
Set D could account for both type-I and type-II models. We 
also notice that Sets A,B have been devised in order to mimic 
the characteristic mass splittings of the MSSM Higgs sector. 



Before coming to grips with the analysis of the Higgs- 
strahlung events in the general 2HDM, it is interesting 
to briefly reconsider the analogue process in the simpler 
framework of the SM. In Fig. [S] we display the total (one- 
loop corrected) cross section (left panel), together with 
the relative radiative correction 8 r (right panel), as a 
function of the SM Higgs boson mass Mh. For complete- 
ness, and for illustration purposes, we have also included 
the scenario corresponding to the last Higgs boson mass 
segment ruled out by LEP at a center-of-mass energy 
\/s = 200 GeV (see the leftmost vertical band in that 
figure). The corresponding production rates for the ILC 
at yfs = 500 GeV and y/s = 1 TcV are smaller than 
in the LEP case due to the suppression of the s-channcl 
amplitude by the Z-boson propagator at higher energies, 
see Eq. ([5]), and the larger mass of the produced Higgs 
boson. Still, the cross-sections for producing SM Higgs 
bosons of a few hundred GeV at the startup ILC energy 
(y/s = 500 GeV) lie at the 40 - 60 fb level, which leads 
roughly to ~ 25,000 Higgs events for an expected in- 
tegrated luminosity of 500 fb -1 . In turn, the one- loop 
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M H [GeV] M H [GeV] 



FIG. 5: Total cross section a(e + e~ -> HZ ) (in fb) at one-loop (left panel) and relative one-loop correction S r (right panel) for 
= 500 GeV and 1 TeV, as a function of Mh and within the SM. In the leftmost side band we have computed also the curve 
corresponding to LEP 200, which falls abruptly from rather high values down till reaching the end of phase space. The thin 
vertical band corresponds to the narrow exclusion region 160 < Mh < 170 GeV determined by the Tevatron collaborations [HEI ]. 



radiative corrections may be either positive or negative, 
depending on the center-of-mass energy, and lie generally 
at the level of a few percent, as can be seen in the right 
panel of Fig. [5] In this panel, we present the evolution 
of the correction parameter 5 r defined in Eq. (|2"Tj). as a 
function of Mh . The three peaks (also barely seen in the 
left panel) are correlated to the production thresholds of 
WW~, Z°Z° and tt pairs. 

In Figs. [BK]wc illustrate the fundamental phenomeno- 
logical features associated to the process e + e~ — > h°Z° 
in the scope of the general 2HDM. Fig. [6] summarizes 
the pattern of radiative corrections S r projected onto the 
(tan/3, A5) and (sin a, tan j3) planes, the former at fixed 
a = f3— 7r/2 and the latter with A5 = —2. We remark that 
for a = (3 — tt/2 the h°Z°Z tree-level coupling takes on 
the SM form and therefore is maximal. For definitcness, 
these plots have been generated for a Higgs boson mass 
spectrum as in Set B (cf. Table UJ, and at the fiducial 
ILC startup center-of-mass energy, yfs = 500 GeV. Al- 
though the range tan/? > 1 is usually the preferred one 
from the theoretical point of view, we entertain the possi- 
bility that tan (3 can be slightly below 1 in order to better 
assess the behavior around this value. Fig. [S] illustrates 
that the allowed region in the (tan /3, A5) plane is severely 
restrained by the theoretical constraints stemming from 
the perturbative unitarity and vacuum stability. On the 
one hand, A5 > values are strongly disfavored by the 
vacuum stability condition; on the other hand, the uni- 
tarity constraints tend to disfavor moderate and large 
values of tan (3 (specially tan j3 values significantly larger 
than 1) as well as of tan/3 <C 1. Altogether these con- 
straints set an approximate lower bound of A5 ~ —10 
for tan/3 ~ 1 and a rigid upper bound excluding almost 



all positive values of A5 for any tan/3. The combined 
set of constraints builds up a characteristic physical do- 
main, with a valley-shaped area centered at tan/3 ~ 1 
that sinks into the A5 < region and becomes narrower 
with growing |As|. 

The fact that the curves of constant 5 r do not depend 
on tan /3 (left panel of the Fig. [5]) is related to the choice 
a = /3— 7r/2, which we have made in order to consider the 
situation where the tree-level cross-section for h°Z° pro- 
duction is maximal. For this choice of the mixing angle 
a, it turns out that no particular enhancement shows up 
for any value of tan /3, neither from the 3H self-couplings 
nor from the Higgs-top quark Yukawa couplings. Indeed, 
under the condition a = j3 — tt/2, the h°tt Yukawa cou- 
pling does not depend on tan /3: 



Vtt 



e cos a m t 
2 Mw sw sin /3 



em t 



a=0- 



2M w s 



IV 



(22) 



and at the same time the 3H self-couplings which are 
relevant for this channel become independent of tan /3; 
notice, for example, that 



Ah°H°H 



2M w s 



w sw 



(Mj 



2M, 



H» ) 



4A 5 M^ s w 



(23) 



The final outcome is that the potential dependence of 
the computed obscrvables on tan /3 vanishes as long as 
we stick to these a = (3 — tt/2 configurations - in which 



10 



.-3 




tan P 

sin a 

FIG. 6: Contour lines for the ratio S r in the [tan/3, As] plane assuming ot = /3 — 7r/2 (left panel); similarly, in the [sin a, tan/3] 
plane for As = —2 (right panel). The dark gray shaded areas stand for the regions excluded by the vacuum stability bounds, 
whereas the light gray shaded areas signal the domains excluded by the perturbative unitarity bounds. The results are obtained 
for y/s — 500 GeV and the Higgs boson masses as in Set B of Table [T] 



the tree-level coupling h°Z°Z is maximum and formally 
equivalent to that of the SM 3 . As a result, the only feasi- 
ble mechanism able to significantly enhance the quantum 
effects in these scenarios is by increasing the value of the 
IA5I parameter (towards more negative values, so as to 
be consistent with vacuum stability). 

Once we depart from the a — /3 — tt/2 setting, we 
recover of course the expected dependence of the com- 
puted observables with tan /3 (cf. right panel of Fig. [5]) , 
but then the lowest order h°Z° production cross section 
becomes smaller. For tan/3 < 1, radiative corrections are 
boosted as a result of the enhanced 3H self-couplings, 
and partially also due to the Higgs-top quark Yukawa 
couplings. Either way, their overall effect is to suppress 
the tree-level signal, as such leading corrections are trig- 
gered primordially by the finite Higgs boson WF cor- 
rections. Another source of enhancement of 8 r appears 
near the region where /3 ~ a. However, this effect is not 
caused by a real increment of the one-loop cross-section 
0-(o+i) ^ j-,ut by a mere suppression of the cross-section at 
the Born-level, i.e. cr' ' in the denominator of Eq. (pH]) . 
and in this sense it is an uninteresting situation. 

In Fig. [7] wc explore the evolution of the cross section 
as a function of the center-of-mass energy. Wc include, 
in each plot, the tree-level contribution and the cor- 
responding loop-corrected value, a^ 0+1 ^ , for different val- 
ues of A5. The right panel tracks the related behavior 
of the quantum correction S r as a function of y/s, for 
the same set of A5 values. The plots are generated for 
Set B of Higgs boson masses, at fixed tan/3 = 1 and 
a = /3 — it/2. The leading-order cross section curve 



3 An equivalent discussion would hold for the H°Z° channel, under 
the complementary condition a = /3. 



exhibits the expected behavior with y/s, as it scales with 
the s-channel Z°-boson propagator, namely proportional 
to l/(s — M§). A similar pattern is also followed by the 
full loop corrected cross-section. 

The range where the relative one-loop correction S r is 
positive is very reduced and it is confined to a regime 
where the center-of-mass energy is around the startup 
value for the ILC (y/s > 500 GeV) and where A5 adopts 
the (small) positive values allowed by the constraints. It 
should not come as a surprise that this behavior is sim- 
ilar to the SM result found in Fig. [5j indeed, being A5 
small there can not be significant 2HDM enhancements 
(not even from tan/3, which is 1) with respect to the cor- 
responding SM process. Let us recall that for Set B the 
maximum value of A5 allowed by the vacuum stability 
bounds is A5 ~ 0.65. Large negative A5 values give rise 
to significant (negative) quantum effects, and translate 
into loop-corrected cross sections er(° +1 ) depleted 40% to 
70% with respect to the tree-level predictions in the en- 
tire range from y/s = 500 GeV to y/s = 3 TeV. Since 
the tree-level h°Z°Z° coupling is equivalent to the SM 
coupling HZ°Z° for a = /3 — tt/2, the right panel of 
Fig. [7J also illustrates the departure of the 2HDM loop- 
corrected cross section with respect to the tree-level SM 
cross-section, namely [o^hdm ~ a SM)l a SM- 

The behavior of cr(e+c- -> h°Z°) MS c.l function of the 
Higgs boson mass M h o is presented in Fig. |5J This fig- 
ure is the 2HDM counterpart of Fig. [5] corresponding 
to the SM case. We superimpose the tree-level and the 
loop-corrected cross-sections for the Set B of Higgs boson 
masses, by setting tan/3 = 1, a = (3 — tt/2 and using three 
different values of A5. The center-of-mass energy is set- 
tled at the fiducial value y/s = 500 GeV. Obviously, the 
raise of the Higgs boson mass M h o implies a reduction of 
the available phase space, so that the cross-section falls 
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FIG. 7: Total cross section a(e + e~ — > h°Z°) (in fb, left panel) and relative one-loop correction 8 r (in %, right panel) as a 
function of yfs using Set B of Higgs boson masses (cf. Table [T]); for tan/3 = 1, a — ft — n/2 and three different values of A5. 
Shown are also (cf. the right vertical axis of the left panel) the number of events per 500 fb -1 of integrated luminosity. 



down. In the left panel of Fig. both the tree-level and 
the loop-corrected cross-sections decrease monotonically 
with the growing of M h o. The Z°Z° and W + W~ thresh- 
olds are also barely visible therein. In the right panel of 
the same figure we observe a steady increase of the nega- 
tive value of the correction for heavier Higgs boson mass, 
whereas the case with positive correction remains almost 
stable. 

Finally, in Fig. [9] we display the tree-level contribu- 
tion and the total (one-loop corrected) cross section (left 
panel), together with the relative radiative correction S r 
(right panel), as a function of the parameter A5. These 
plots have been generated in the same benchmark con- 
ditions as in Fig. [S] and using a fiducial ILC start-up 
center-of-mass energy s/s = 500 GeV. Let us remem- 
ber that the red shaded area on the right hand side of 
that figure (specifically for A5 > 2.4) stands for the re- 
gion excluded by the vacuum stability bounds, whereas 
the light gray shaded area on the left (As < —10.5) sig- 
nals the domain excluded by the perturbative unitarity 
bounds. The cross section is seen to follow the expected 
behavior, which we have derived from the estimate of the 
leading effect - cf. Eq. (fl"6|) - in combination with the 
structure of the trilinear coupling (|23[) , namely a negative 
quadratic dependence of the scattering amplitude, which 
transates into a(e + e~ — > h°Z ) ~ (a — &A|) 2 , emerging 
ultimately from the finite WF corrections to the external 
Higgs boson leg (notice the inclusion of the leading quar- 
tic corrections). For moderate negative A5, the relative 
size of the one-loop quantum corrections can reach up to 
~ —50%. Positive 5 r , however, can only take place for 
A5 close to zero. 

A rather comprehensive survey of the predicted cross- 
sections for different values of the tree-level coupling 
and different Higgs mass setups is presented in Ta- 
ble |TTJ We display the results for the one-loop corrected 



cross-sections (in fb), together with the relative radia- 
tive corrections (cf. Eq. (|2~Tj) ). for both neutral Higgs- 
strahlung channels h°Z° and H°Z°. Once more, we set 
y/s = 500 GeV; and we work at tan /3 = 1 and maximum 
allowed | Ag [ , since we are mostly interested in spotlight- 
ing the imprints of the 3H self-couplings in the quantum 
effects associated to the Higgs-strahlung mechanism. Let 
us also recall in passing that a = /3 — tt/2 maximizes the 
tree-level h°Z°Z° coupling, whilst H°Z°Z° is optimal in 
the complementary regime, a = {3. From the table we 
may read out that radiative corrections in these regimes 
are certainly large, regardless of the details of the cho- 
sen mass spectrum, the tree-level coupling and the actual 
channel under consideration. In a nutshell: the charac- 
teristic signature of the enhanced 3H self-couplings in the 
pattern of quantum effects on the Higgs-strahlung pro- 
cesses is rather universal and manifests in the form of a 
substantial depiction of the tree-level signal (typically in 
the range of 6 r ~ —20/ — 60%). Fortunately, even un- 
der such a dramatic suppression the final loop-corrected 
cross sections stay at the level of a few tens of fb - thereby 
amounting to a non-negligible yield of ~ 10 3 — 10 4 events 
per 5001b- 1 . 

A note of caution should be given at this point in re- 
gard to some particular configurations of the H°Z° chan- 
nel. It turns out that certain choices of Higgs boson 
masses (e.g. for MSSM-like mass splittings) may lie close, 
or simply beyond, the kinematical threshold for the de- 
cay process H° — > h°h°; that is to say, M H o ~ 2M^o. 
If we then consider a regime wherein a = /3 and sizable 
I A5 1 , then the H° WF-correction undergoes a remark- 
able boost, which results from the combination of two 
independent sources of enhancement, namely i) the ac- 
tual strength of the H°h°h° coupling - which is max- 
imally enlarged in this regime; and ii) the kinemati- 
cal enhancement due to the vicinity of the H° — > h°h° 
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FIG. 8: Total cross section cr(e + e — > h Z ) (in fb, left panel) and relative one-loop correction 5 r (right panel) at fixed yfs = 500 
GeV as a function of M h o for Set B of Higgs boson masses, cf. Table [I] Shown are the results obtained within three different 
values of As and for a = /3 — tv/2. 
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FIG. 9: Total cross section cr(e + e~ — > h°Z°) (in fb) at tree-level and at one-loop (left panel), together with the relative one-loop 
correction 8 r (right panel, in %) as a function of A5. The results are obtained for tan/3 = 1, a = ft — n/2, ^/s — 500 GeV and 
the Higgs boson masses as in Set B of Table [I] 



threshold, which is reflected as a sharp peak in 8 Z a o ~ 
B' (M^ ,M^ ,Ml ). However, these H° WF-corrections 

become so large that the perturbative formula (fT6]) is 
no longer valid and we have to keep the original (unex- 



panded) expression Z^ 2 = 1 + $leY, / H o(M^ i0 ) 

of such corner in the parameter space, the higher-order 
effects involved in the above rcsummed form of the WF- 
corrections become harmless and totally inconspicuous 
and can be safely discarded. If, alternatively, we were 
considering a mass setup such that A/ H o > 2M h o, then 
the decay H° — > h°h° would be open and, in the scenario 
of maximum H°h°h° coupling, it would furnish a very 
large width for H°. In this particular setup, the process 



-1/2 



Out 



would effectively boil down to the double Higgs-strahlung 
channels e + e~ — > H°Z° — ^ h°h°Z° previously explored in 
the literature f22j . 

The following comment is in order to clarify the role 
played by the remaining (potential) sources of enhance- 
ment. We have seen that the various constraints restrict 
very significantly the range of allowed values of tan /?, 
and enforce it to stay around 1 when |As| is maximum in 
the region A5 < 0. As a result, the contributions from the 
top quark and bottom quark Yukawa couplings cannot be 
augmented in the domain where the trilinear couplings 
are maximal. Remarkably enough, let us also emphasize 
that they cannot be enhanced even in the region where 
I A5 1 is small or zero. To see why, notice that although in 



13 





h u Z u 


H U Z U 


a = f3-n/2 


a = /3- tt/3 


a = p — 7r/6 


a = P -7r/3 


a = P — 7r/6 


a = P 


Set A 


(Jm n -v ID 

mux L J 


42.61 


27.89 


8.17 


5.36 


17.86 


10.70 


S r [%\ 


-22.79 


-32.62 


-40.76 


-40.92 


-34.40 


-70.51 


Set B 


Umax [fb] 


24.80 


18.84 


7.24 


7.05 


15.51 


14.99 


Sr [%J 


-53.18 


-52.57 


-45.29 


-41.19 


-56.87 


-68.74 


Set C 


Umax [ft>] 


28.81 


21.97 


7.75 


7.83 


20.82 


25.60 


5 r \%\ 


-43.91 


-42.97 


-39.66 


-36.96 


-44.13 


-48.48 


Set D 


Umax [ft>] 


40.77 


29.25 


9.30 


8.56 


23.95 


28.59 


8 r [%\ 


-17.95 


-21.52 


-25.15 


-15.14 


-20.84 


-29.12 



TABLE II: Maximum total cross section cr (0+1) (e + e~ -> hZ°), for h = h°,H°, at yfs = 500 GeV, together with the relative size 
of the radiative corrections 5 r , for the different sets of Higgs bosons masses quoted in Table [I] The results are obtained at fixed 
tan/3 = 1 and different values of a, with As at its largest negative attainable value, namely: As ~ —9 (for Set A), As ~ — 10 
(Sets B,C) and tan/3 ~ -8 (Set D). 



such region the parameter tan /3 can be very large or small 
and still be compatible with the constraints (cf. Fig. [6]), 
then all the terms in the trilinear couplings (cf. Table 
II of Ref. [32| ) which are not proportional to A5 become 
of order one (i.e. cannot be promoted to high values) in 
the regime where the tree-level production cross-section 
for h° or H° is optimal. Not only so, in the very same 
regime the Yukawa couplings of both the top and bot- 
tom quarks cannot be enhanced either. One can easily 
check all these features explicitly by observing that, in 
the regions where the corresponding tree-level processes 
are maximal, the tan /3-cnhancements arc canceled in all 
the couplings. We have checked numerically that the 
residual corrections (positive and negative) are merely of 
a few percent for any value of tan/3. At the end of the 
day, we conclude that the only sizeable and eventually 
measurable quantum effects on the processes under study 
are those stemming potentially from the enhancement of 
the A5 coupling, not from the Yukawa couplings. 

Let us finally emphasize that the complementarity be- 
tween the neutral final states h Z°/h°A and H°Z /H A°, 
i.e. processes (JTJ and ([2]), is a unique chance for analyzing 
potentially big correlations between quantum effects in 
the 2HDM production cross sections, as there is no sim- 
ilar opportunity in the charged sector. Indeed, the final 
state charged counterparts H ± W T are suppressed owing 
to the fact that the ZH ± W =F vertices are forbidden at 
the tree-level in any 2HDM extension of the SM [56| and 
therefore they can only be studied in more general ex- 
tensions of the Higgs sector WR. or through loop- induced 
H±W T vertices in the 2HDm|1] and in the MSSM [H]. 
In these loop-induced mechanisms, the cross section for 
the associated H ± W =F production in a linear collider is 
rather meager - generally below 1 fb at the startup ILC 
energy (for charged Higgs masses comparable to the ones 
we have considered for the h°Z /h°A production) and 
within the region allowed by the current constraints. We 
therefore deem quite difficult to use the H ± W T chan- 
nel to extract additional information. In our opinion, 
the main task should be focused on performing preci- 



sion tests of the h Z /h°A cl final states, together with 
the H Z°/H°A ones (if the mass of the heavy CP-even 
Higgs is not too large). 



IV. DISCUSSION AND CONCLUSIONS 

In this article, we have concentrated on the analy- 
sis of the production of neutral Higgs bosons in asso- 
ciation with the Z° gauge boson at the future linac 
facilities within the framework of the general (non- 
supcrsymmctric) Two-Higgs-Doublct-Model (2HDM) . 
Our basic endeavor has been twofold: i) on the one hand 
to study the impact of radiative corrections to the final 
predicted rates; and ii) on the other hand to correlate 
such quantum effects with the enhancement potential of 
the 3H self-interactions, which are a genuine dynamical 
feature of the 2HDM - in the sense that it is unmatched 
to its supcrsymmctric counterpart. The upshot of our 
analysis singles out sizable (although negative) quantum 
effects which are correlated to the enhancement prop- 
erties of the 3H self-interactions in the general 2HDM. 
Such large effects are identified fundamentally in the re- 
gion of parameter space with tan j3 ~ 1 and where the 
coupling IA5I is at its maximum possible value compat- 
ible with the various constraints. The quantum effects 
reach typically 8a /a ~ -20%/ - 60%, for |A 5 | - 8/ - 10 
(A5 < 0), this being the crucial parameter that tunes the 
actual size of the 3H self-couplings. Let us stress that the 
most stringent limits on A5 are placed by the conditions 
of perturbative unitarity and vacuum stability 4 . 



4 Recently, a combined analysis of different B-meson physics con- 
straints over the 2HDM parameter space suggests that values of 
tan P ~ 1 could be disfavored for charged Higgs masses Af H ± too 
near to the lowest mass limit of 300 GeV (and certainly below) 
[iioll . However, the level of significance is not high and, in ad- 
dition, our leading quantum corrections are basically insensitive 
to the charged Higgs mass. Therefore, a shift of M H ± slightly 
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Although the vertex corrections to the h°Z Z /H Z°Z 
couplings are sensitive to such 3H self-interactions 
through Higgs-mcdiatcd one-loop diagrams, the most sig- 
nificant quantum effects are induced by the wave-function 
renormalization corrections to the h°/H° Higgs boson ex- 
ternal lines. These (negative) effects triggered by the tri- 
linear couplings are of order a cw \^ H , and are therefore 
very responsive to changes in the value of the param- 
eter A5 in the Higgs potential of the 2HDM. In turn, 
the gauge-boson and fermion one-loop contributions re- 
main subleading (including the effects from the Yukawa 
couplings) as they show no remarkable departures from 
their SM analogues in the relevant regions of parame- 
ter space - which yield 5a /a at the level of a few per- 
cent. Most important is also the fact that the combined 
analysis of such Higgs-strahlung events with the previ- 
ously considered Higgs-pair production processes (h°A°, 
H°A°) - see [IH for details - could be instrumental as a 
highly sensitive probe of the underlying architecture of 
the Higgs sector. At the fiducial center-of-mass energy 
value of y/s = 500 GeV, and when the genuine enhance- 
ment mechanisms of the 2HDM are active, the hAP events 
are remarkably strengthened at the one-loop level while 
the hZ° channels arc simultaneously depleted. 

We can assert that the described pattern of leading 
quantum effects emerges as a kind of "universal" fea- 
ture of the 2HDM as far as the predictions for the 
Higgs-strahlung processes are concerned, meaning that 
this pattern is virtually independent of the details of 
the Higgs mass spectrum and of the actual production 
channel (h°Z°, H°Z°). The same is true for the pair- 
wise production channels [HI . Therefore they all can be 
physically significant provided the produced Higgs bo- 
son is not too heavy for the actual center-of-mass energy. 
Focusing on the Higgs-strahlung processes, the rise of 
large (and negative) radiative corrections to their cross- 
sections within the 2HDM can be regarded as a char- 
acteristic imprint of a general two-Higgs-doublet model 
structure. Furthermore, the presence of a Z°-boson in 
the final state (and hence of its clear-cut leptonic sig- 
nature Z° — > l + l~) should enable a rather comfortable 
tagging of these Higgs events following a method entirely 
similar to the original Bjorken mechanism (35| in the SM. 
Although the corresponding cross-sections at the highcr 
operating energies planned for the future linear colliders 
are smaller (as compared to LEP), they are neverthe- 
less sufficiently sizeable (typically in the range 10 — 40 
fb) for a relatively comfortable practical measurement 
at the higher luminosities scheduled for these machines. 
For example, the associated decay of the accompanying 
neutral CP-even Higgs-boson would manifest basically in 
the form of either i) collimated and highly-energetic b- 
quark or r-lepton jets, for M h o < 2 My ~ 180 GeV; 



upwards restores the possibility of tan ~ 1 at, say, la without 
altering significantly our results (as we have explicitly checked). 



and ii) h°,H° — > W + W — > l + l + missing energy; or 
H° -> Z°Z° -> l+l'l'+l'-, for M H o > 2 My (in the opti- 
mal regime for the tree- level process). 

Some discussion on how the corresponding MSSM 
Higgs boson production cross-sections compare to the 
2HDM ones seems also appropriate. In the MSSM, the 
maximum quantum effects on hZ° production are typi- 
cally milder. Here, the dominance of the WF-corrections 
is also the main source of one-loop effects [29[ . They are 
usually reabsorbed into the tree- level ZZ/i couplings ((H), 
specifically in the CP-even mixing angle a, which then 
becomes an effective mixing parameter. Consider, for 
example, the behavior of the characteristic trigonome- 
tric coupling cos 2 (/3 — a) in the ZZH° interaction vertex, 
which in the MSSM dies away with growing values of 
Mpa (both at the tree- level and at one- loop). It follows 
that in a situation where Mpp is sufficiently heavy, say 
above 200 GeV, the MSSM would predict measurable 
rates for the channel h°Z° (and also for H°A°, if y/s is 
sufficiently high), whereas the complementary channels 
(H n Z°/h A ) would be virtually below the observability 
threshold. In contrast, in this particular SM-like regime 
(in which sin(/3 — a) —> 1) the h°Z° cross-section in the 
2HDM could well be exhibiting the trademark suppres- 
sion induced at one-loop by large 3H self-couplings. 

Furthermore, in the MSSM, the dynamical origin of 
the leading WF-effects does not reside in the Higgs bo- 
son self-couplings, but on the Yukawa couplings with 
fcrmions and also in the large Yukawa-like couplings of 
Higgs bosons with sfcrmions, particularly with the stop. 
Thus, in both models (2HDM or MSSM) the main source 
of quantum effects emanates from the renormalization of 
the Higgs boson external lines, but the kind of interac- 
tions involved in each case is radically different. On the 
quantitative side, the impact of the quantum effects in 
the MSSM case turns out to be considerably milder and 
as a consequence the cross-sections for Higgs-strahlung 
production remain similar to the SM cross-section (i.e. 
for the Bjorken process). The corresponding 2HDM 
cross-sections, instead, can be significantly smaller - if 
taking the same Higgs boson masses (e.g. Sets A and B 
of Table HI which mimic the MSSM Higgs boson mass 
spectrum) owing to the aforementioned large suppres- 
sion effect from the Higgs boson WF renormalization. At 
the same time, as we have already mentioned, the cross- 
sections for neutral Higgs boson pair production in the 
general 2HDM - cf. Eq. ([1]) - become substantially larger 
than the MSSM counterparts (for the same Higgs bo- 
son mass spectrum), carrying very sizeable and positive 
quantum effects [32J. Therefore, an interesting combined 
picture emerges in the 2HDM context, in which a large 
(~ 50%) enhancement of the pairwise neutral Higgs bo- 
son production is simultaneously accompanied by a dras- 
tic drop (of similar size) of the Higgs-strahlung events. As 
this situation is completely impossible to realize in the 
MSSM, this feature could be used as a strong character- 
istic signature to discriminate between these processes in 
the general 2HDM and in the MSSM. In other words, it 
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could provide an essential quantum handle enabling us 
to perform a proper identification of the kind of Higgs 
bosons produced in a linear collider. 

In summary, we have analyzed the classical Higgs- 
boson strahlung processes at linear colliders in the light 
of the general 2HDM, and we have elucidated very sig- 
nificant quantum imprints which spotlight, once more, 
the stupendous phenomcnological possibilities that triple 
Higgs boson self-interactions could encapsulate in non- 
supersymmetric extended Higgs sectors. To be sure, re- 
gardless of whether the LHC is finally capable to dis- 
cover the Higgs boson(s), a paramount effort will still 
be mandatory in order to completely settle its experi- 
mental basis and to uncover the nature of the spinless 
constituents behind the electroweak symmetry-breaking 
mechanism. As we have shown, experiments at the fu- 



ture linear collider facilities could play a crucial role in 
this momentous endeavor. 
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